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Abstract 

The authors propose a robust semi-parametric empirical hkehhood method 
to integrate all available information from multiple samples with a common 
center of measurements. Two different sets of estimating equations are used 
to improve the classical likelihood inference on the measurement center. The 
proposed method does not require the knowledge of the functional forms of 
the probability density functions of related populations. The advantages of 
the proposed method were demonstrated through the extensive simulation 
studies by comparing mean squared error, coverage probabilities and average 
length of confidence intervals with those from the classical likelihood method. 
Simulation results suggest that our approach provides more informative and 
efficient inference than the conventional maximum likelihood estimator when 
certain structural relationships exist among the parameters for these relevant 
samples. 
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1. Introduction 

A common problem in clinical trials and medical research is how to accu- 
rately and efficiently estimate parameters of interest when the current sample 
size is small due to cost and time constraints. Usually there might exist cer- 
tain surrogate populations with low sampling cost that could provide relevant 
information for the population of direct inferential interest. In this article, 
we propose a robust semi-parametric method to integrate related information 
from different sources to improve the classical likelihood method. 



Preprint submitted to Elsevier 



October 3, 2012 



The classical likelihood approach is arguably the most widely used method 
in statistical inference. It has been routinely applied in almost all the statis- 
tical applications. Despite the great success and excellent asymptotic proper- 
ties, the classical likelihood has limitations associated with making inference 
for small sample sizes. Consider a thought experiment as follows. Suppose 
that a random experiment is to toss a coin twice. The parameter of inter- 
est, denoted as 6i, is the probability of turning up head for this coin. The 
maximum likelihood estimator (MLE) of 6i is denoted as 6i. If the coin is 
a fair one, the MLE will obtain the following P{6i = or 1) = 1/2 and 
P{6i = 1/2) = 1/2. Thus, one would have 50% chance to make a nonsensical 
decision by using the MLE when the sample size is only two. In addition, 
suppose that for some reason we cannot use this coin any more but we can 
flip another coin instead. In this situation, the classical likelihood approach 
would not consider the second experiment since it comes from a different 
population unless a functional relationship between the two parameters is 
known. If the second population is related to the first one due to some un- 
known link between these two parameters, one should be able to utilize this 
connection and make better statistical inference. 

Different statistical methodologies have been proposed in the literature to 
integrate information from different sources (or populations) in a very general 
setting, see Wang, van Eeden and Zidek (2004), Fu, Wang and Wu (2009) 
and referees therein. Most of these methods, however, face the challenge of 
accurately validating or evaluating the relevance of all related information 
to handle the possibility of introducing a significant bias or contaminating 
the current sample. In other words, the magnitude of integration must be 
controlled carefully and in addition likelihood weights must be chosen judi- 
ciously in order to achieve any desired improvement in statistical inference. 
We propose to tackle this difficult problem using a robust semi-parametric 
empirical likelihood method to gain more accurate and robust inferential 
results. 

Empirical likelihood, which was first introduced by Owen (1988), is a 
nonparametric method of inference based on a data-driven likelihood ratio 
function. It allows the statistician to employ likelihood methods, without 
specifying a parametric model for the data. It enjoys both the fiexibility 
of nonparametric methods and the efficiency of parametric likelihood. As 
shown in Qin and Lawless (1994), empirical likelihood is a prominent effi- 
cient tool in estimating parameters by incorporating estimating equations 
into constrained maximization of the empirical likelihood function. In the 
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problem we consider, the relevant information from different sources could 
be used by incorporating extra set of estimating equations in the empirical 
likelihood framework. 

To obtain robust estimates, we use median as an estimate of center instead 
of mean. We propose here using two different kinds of estimating equations; 
one uses median and the other one uses a smoothed version of median. The 
smoothing technique is the one proposed by Shi and Lau (1999) to improve 
the coverage accuracy. Our method can be easily generalized to multiple 
samples with relevant information. Without loss of generality, we consider 
data with two populations. 

The rest of the paper is organized as follows. The methodology frame- 
work, the proposed empirical likelihood approach, and its theoretical proper- 
ties are presented in Section |2l Results of simulation studies demonstrating 
the empirical performance are provided in Section 3. Conclusion and some 
discussion are provided in Section HI 

2. Methodology 

Suppose there are two groups of data from different population but shar- 
ing the same parameter of interest. Assume that 

The second group of data y^, . . . , y^^ might be different from the first pop- 
ulation, and 

yi,- -,yn2 ~ h{y,^)- 

Our goal is to estimate by using both samples. Directly using the log- 
likelihood 

ni n2 

^log/(x,,,0) + ^logMy„0), 

1=1 i=l 

we might get a biased estimation due to the difference between the two 
populations. 

We propose a semi-parametric empirical likelihood method which only 
requires the independence of these two samples. To combine the second 
sample with the first one, we use the following semi-parametric empirical 
likelihood 

rii n2 
^ = ^^ogf{xi,e) + ^logn2Pi, 

i=l i=l 



3 



where 

n2 n2 

Pi > 0, = 1, "^PiQiVi, 0) = 0, 

i=l i=l 

and giVi, 9) is an estimating function. From the empirical hkehhood theory, 
we know Pi is maximized by l/{n2 [l + X^g{y^, 0)] }, where A is the Lagrange 
multipher. We can rewrite the log likelihood function as 

ni n,2 

m = 5^ log f{x,, e)-J2 log [1 + ^^9{y^, 0)] , (1) 

i=l i=l 

and 

= argmax£(0). 


We call the robust semi-parametric empirical likelihood estimate (RSPELE). 

The advantage of the log profile likelihood function is that it does not 
depend on the likelihood weights which could be difficult to choose. In our 
propose method, we do not require that the probability density function 
of the second population is identical to the first population. By using the 
empirical likelihood method, we do not even need to specify the functional 
form of the underlying distribution of the second population. Therefore, we 
can gain robust estimates in the sense that model mis-specification problem 
is avoided. Consequently, our method can be employed in a relatively wide 
range of applications when the functional form of the probability density 
function is not known. 

In the following, the theoretical properties of the proposed RSPELE es- 
timator will be presented. For clarity, all proofs are postponed to the Ap- 
pendix. Theorem [T] below shows that under some regularity conditions, the 
RSPELE estimator is consistent to ^o- 

Theorem 1. Let Xi, . . . , a3„^ be i.i.d. from f{x, 0) and tfi, . . . , i-i-d. 

from an unknown distribution. Assume that f{x,0) satisfies the regularity 
conditions given in Shao (2003) on the normality of the maximum likelihood 
estimator in parametric models. Let 0o be the true parameter. We further 
assume that 

(Al) There exists a matrix > such that 

E[g{y^^,0)] =^(0-0,) + o{\\0-0o\\) 

as - 00 ^ 0. 
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(A2) ({6) = E[\\g{yi,0) - g{y^,eo)\\'^] exists when \\e - OqW is sufficiently 
small and is continuous at 6q 

(A3) E[giy,,eo)g^{y„eo)] > 0. 

Then, it follows that 6q in probability in the neighborhood of 6q such 
that{0: ll^-^oll <n~^/^}. 

The asymptotic distribution of the 6 and A is shown in Theorem [2l 

Theorem 2. In addition to the conditions of Theorem U\, we assume that 

ui dg(y,6) 

— b, where b is a constant. We also assume that — - — exists with 

n ad 
probability one and the set of its discontinuity points has zero probability. 

Then 



V^{0 - Oo) 



iV(0,^i), 
iV(0,52), 



where 
So 
Si 

S2 



bI{do) + (1 - b)E 



dg 
W 



bS^'li0o)S^' + (1 - b)S^'E 



b^-^E 



dg 



s^h{e,)s,'E 



dg 

dg 
80^ 



dg 
80^ 

T 



30' 



T 



-$-1 + (1 - b)^-'E S,'E 



$-1 



dg 



00' 



$-1 



x<l> 
Egg 



+ ( %) S^^E ( % 



d0' 



80' 



T 



$-1 



T 



(2) 



and I{0o) is the Fisher information about 6o contained in X . 

The asymptotic distribution of 2[i{0, A) — i{0o, Aq)] is given in Theorem 

El 



Theorem 3. Assume that the assumptions made in Theorem [H hold. The 
limiting distribution of2i{0,X) — 2i{0Q,XQ) is the same as the distribution 
of 



Ui 
U2 



V 



-1 



^0 

b)E 



'l-b)E 



m 



d9_Y\ 
dO^J 

J 



Ui 

U2 



under Hq : = 6q, where Ui is independent ofU2, 

U,^N{0,bI{6o)), t/2~iv(0, {l-b)E[g{Y,eo)9^iY,e, 

and 



( 



0-1 
'-'0 

dg 



where V^'^'^^ 



._l_<|.-i + $-1^ 



dg 



dg 

w 

y{22) 

dg 
86 



-A 



T 



We note that other test statistics, for example, a test statistic based on 
Theorem may also be used. 

Estimating equations provide a very flexible way to specify how the pa- 
rameters of a statistical model should be estimated. They serve as constraints 
in maximizing the empirical likelihood. Qin and Lawless (1994) showed 
that the empirical likelihood method is an efficient tool for point estima- 
tion through estimating equations. In this section, we consider two different 
kinds of estimating equations using the information of median, since median 
is robust with respect to the outliers, one may use 

g^{y,,0) = l{y,<m{0))-l/2 

as estimating function based on the second group data, where m{6) is the 
median of /(a;, 6) and l(-) is the usual indicator function. It is easy to verify 
thatE(^i(t/„6>))=0. 

Due to the discontinuity of gi, we may use the smoothed version of the 
constraint which was motivated by Shi and Lau (1999). First of all, we define 
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the estimating equation for the smoothed empirical hkehhood. In general, 
let K be the r^'^-order kernel (Shi and Lau, 1999), such that 

r r 1, ifj = 0, 

/ u^K{u)du = I 0, if 1 < J < r - 1. 
J y Cr, if j = r, 

where r is a positive integer. Define ip{u) = j K(t)I{t < u)dt. For any /i > 0, 
let iph{u) = ip{u/h) where h is called the smoothing parameter. The kernel 
K is a symmetric probability density with bounded and compact support. 
Let the estimating function g2{y,9,h) = iph{m{6) — y) — 1/2. Therefore, 
g2{y, 0, h) is continuous with respect to y, but it is not a fixed function as the 
smoothing parameter varies. See Shi and Lau (1999) for details. In addition, 
by using the arguments similar to those stated in the Theorems [T] to [3] and 
Shi and Lau (1999), we may get the similar asymptotic results. 

3. Numerical Experiments 

3.1. Data Fusion with Conventional Empirical Likelihood 

Simulation studies are carried out by performing data fusion when two 

samples are available. The first sample X = (Xi,...,X„J is generated 

from standard normal distribution and the second sample Y = (Yi, . . . , F^j) 

is generated from normal, double exponential, or t-distribution respectively. 

The sample size of first sample, ni, is 10 and for the second sample, the 

sample size, n2, varies from 10, 20 to 30. 

First of all we use the median constraint, so the log likelihood function 

of the simulation model is 

/(/.) = - 5^ log ^2^sl + ^^^-^ - E {1 + ^[^(^^ ^ - 1/2]} , 
i=i ^ 1 ^ i=i 

where sf is the MLE variance, = Y17=ii^i ~ ^)^- 

We present the mean square error (MSE) ratio of RSPELE to MLE based 
on 1,000 replications in Table [H The simulation results show that RSPELE 
performs well except in the situation when t second population is normally 
distributed with large variation as the first one. When the second sample 
size is increasing, the RSPELE becomes more accurate. Moreover, we have 
smaller MSE of RSPELE when the data of the second population is more 
concentrated around the center, for example, the double exponential distri- 
bution. 
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Table 1: MSE Ratio of RSPELE to MLE based on 1,000 replications. 

Distribution of sample 2 
712 A^(0,l)iV(0, 1.25)iV(0, 1.5) A^(0,2)iV(0,3) ^3 DE{0,0.5) DE{0,1) DE{0,1.5) 

10 0.776 0.857 0.929 1.017 1.028 0.787 0.747 0.431 0.717 0.875 

20 0.569 0.693 0.769 0.871 0.956 0.621 0.580 0.209 0.515 0.705 

30 0.453 0.559 0.661 0.789 0.926 0.462 0.467 0.137 0.387 0.582 



3.2. Smoothed Empirical Likelihood 

In this section we demonstrate the smoothed version of the estimating 
equations. The kernel we chose is the same as the one used in Shi and Lau 
(1999), 

[ otherwise. 

The simulation model is identical to the first experiment. Four values of the 
smoothed parameter are used which are n2 to the power of -1, -3/4, -1/2 and 
-1/4. The log likelihood function of the simulation model is 

Kf^) = - E V 27rs? + ^^^^^ - E log {1 + Xg^im, e,h)}, 
i=i ^ ^ ^ i=i 

where si = — X]r=i(-^« ~ •^)^ MLE variance. 

We provFde the MSE ratio of RSPELE to MLE based on 1,000 rephca- 
tions in Table [H Results of the smoothed version are slightly better than 
the results of the median version, no matter which smoothing parameter is 
chosen. When the underlying distribution of the second population is not 
the same as the first population, the RSPELE estimate performs better than 
the MLE. When the sample size of the second population is increasing, the 
RSPELE estimate is more accurate. 

3.3. Confidence Intervals 

In this subsection, we construct the confidence interval for the median 
by bootstrapping. In this simulation study, the first sample X is generated 
from standard normal distribution and the second sample Y is generated 
from normal, double exponential, or t-distribution. The sample size of X is 
10 and of Y varies from 10, 20, to 30. The size of the bootstrapped sample 
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is 200 and the number of iterations is set to be 1,000. First of all, we use the 
median estimating equation and record the coverage probabilities and the 
simulated average confidence interval lengths (AL) in Table H] for nominal 
levels of 80, 90, 95, and 99 percent. The coverage probabilities and the AL of 
using the smoothed version of the estimating equation are recorded in Tables 
Eland [HI with different smoothing parameters which are n2 to the power of -1 
and -1/2. We report the results of MLE in Table O Since the results of MLE 
do not depend on the second population, we further compare the coverage 
probabilities and AL as in Tables H] to |6] with Table [31 

The results of smoothed version are better than median version in terms 
of the coverage probabilities. The coverage probabilities of RSPELE and 
MLE are very close but the confidence intervals of RSPELE are about 10% 
narrower than of MLE. The results of RSPELE when the underlying distri- 
bution of the second population is either t or double exponential distribution 
are better than the results of RSPELE when underlying distribution is nor- 
mal distribution. That is because normal distribution is flatter than t and 
double exponential. Consequently, if the second population provides a good 
information about the center we can use it to get better estimates. 

4. Discussions 

In this paper, we propose a robust semi-parametric empirical likelihood 
in a multiple-sample model with common measurement of center. We use 
two different kind of estimating equations of information about the median. 
Simulation studies have shown that the second population could provide very 
useful information on the parameter of interest by comparing the performance 
of various commonly used measures for evaluations. 

Acknowledgements: We would like to thank Dr. Jing Qin for his 
generous and valuable comments and suggestions for this article. 
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Appendix: Proofs 

Proof of Theorem [I] 

We rewrite tlie equation (HD as i{d) = ii{0) + i2{0), wliere ii{0) 
ZZi ^ogfixi, d) and l2{e) = - J:Zi log[l + >^^9{y^, 0)]. We denote 

N{Oo) be the neighborhood of Oq such that {6 : \\6 - 6o\\ < n^^/^}; 
ON^Oq) be the boundary of N{Oo), i.e. all 0* such that — ^oll = n 
A^int(^o) be the neighborhood of 9o such that {6 : \\6 - 6o\\ < n"^/^}; 
n = ni + n2. 

Til 

Case 1. 0<6i< — < b2 < oo, where bi and 62 are two constants. 
n2 

In view of the proof of Theorem 4.17 of Shao (2003), it follows that for 
any e > 0, 



-1/2. 



P 



(^o) < for all e 9iV(6>c 



> 1-^, 



(3) 



for large ni. 

By Assumptions (Al)-(A3), applying a similar approach as in Owen 
(2001), it can be shown that A = 0(n^/^). In light of Bai, Rao and Wu 
(1992), Qin and Lawless (1994), and under the assumptions (Al)-(A3), it 
follows that for any e > 0, 



{6. 



?2(6'o) < for all 6», G 9A(0o) 



> 1 



(4) 



for large n2- 

Combining (|3]) with 



we have for any e > 0, 
(6>,) - £(0o) < for all 6>, G ON^e^)] > 1 - £, 
for large n. Therefore, there exists G Aint(0o) such that 



80 



and = arg max ^(^)- 



By the definition of N{0q), it follows that 6 ^ 6q in probability, 
ni n-i 
Case 2. )■ 00 or )> 0. 

n2 n2 

The consistency of 6 can be shown similarly. The details are omitted. 

□ 
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Proof of Theorem [H 
We denote 

dh{e) ^d\ogf{x,,e)^ 



d9 89 ' 



A; 



Since 

^^ = and g2„^(^,A) = 0, 
by applying Taylor's expansion, it follows that 
= w^,{e) + q,^^i0,X) 

^ni[0o) H T^Tf - ^O) H TTT + Op{On), 

oO oX 



and 

= q2n,{0,X) 



2n2lt^0, H {i> - ifQ) H A + Op[ 



where (5„ = ||0 — 0o|| + 11 All. It is noted that 









0) 


dX^ 






0) 


do"" 






0) 



E 



dOde 

1=1 



T ' 



n2 



1=1 

dg{yi,0o) _ 

T 1 

uu 

^^2 



n2 

^ dO 

1=1 



dX 



T 



i=l 
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Hence we have 

f _dWnM gqin,(^0.0) \ 

which can be written as 



0-00 
A 



Wnji0o) + Op{dn) 
(00,0) +0p{6r. 



V^{0 - 0o] 



=^„,(0o) + oJ-= 



1 f 8 



where 



/ l dwn,{0o) l gqi„,(6>o,0) 



It can be shown that 

5n = Op(n-l/2); 

n-^I^WnA0o) ^ N{Q,hI{0o))\ 

-n-'/'q^n, (00, 0)^N (o, (1 - b)E [g{Y, 0o)g^{Y, 0o)] ) ; 
1 dwn,{0o) a.s. 



n 80^^ 

ldqirJ0o^ a.s. 



bl{0oy, 



l dq,^^i0o,O) a^ _ f dg{Y,0o) \ 
n 80^ ^ ' \ 80^ )' 



l 8q,^^{0o,O) 
n 8X^ 
Therefore, it follows that 



il-b)E{g{Y,0o)g^{Y,0o)). 



^{0 - 6>o) 
y/nX 

where W is the covariance matrix of 

12 



iV(0,H^), 



with 



V 



( 



{\-h)E 



dg{Y,eo) 



{l-h)E 



il-b)E[g{Y,eo)g^{Y,eo))] ^ 



^21 V22 



Ui is independent of U2 and 

U,^N{0,bI{eo)), C/2~iv(o, {l-b)E[g{Y,0o)g^{Y,6o)] 

Since 



"^11^^12^2.1 



^22.1 



/y(ll) y{12)\ 
[Yi21) Y(22) j 



"^2^^21^11.2 ^22^ + ^2"'^^21^11.2^12V'22^ 

with V11.2 = Vii — V^i2V'22^V2i and V22.1 = V22 — V2il^n^^i2; it follows that 
where 

Proof of Theorem 

Assume that the assumptions made in Theorem [3] hold. A statistic for 
testing Hq : 6 = 6q is given by 

{rii n2 
log fix,, o)-Y^ log [1 + ><^giy^, 0) 
1=1 i=l 
( n2 112 

-2\j2^ogf{x,,eo)-J2^og\l + X^g{y,,do) \. (5) 



i=l 



i=l 
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Expending £(0, A) around (^o; 0) by Taylor's expansion, we have 

rii 712 



i=l 
ni 



1=1 



Viog/(a.„ Oo) + {6- eofw^M + ke - Oof ^'":'}^'\ e - e. 



1=1 



-AV(«..0)-iA^*"^"'*'°''" 



2- ^ + 



n 



+ Op(l). 



92n2(^0,0) 



V^-Miag(V^ll ^22)^^"' 



n 



l^-^^2„2(^0,0)^ 



(6) 



Under Hn, we have 



"2 

-E 



1 + Ao 6>o 



0, 



which imphes that 

H "2 



«=1 



^ "2 



-(l-6)V22'-q2n2(^0,0)+O,(l). 

^^2 
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Hence, 



2 log [l + A^g(y„0o 

i=l 



i=l 



i=l 



^2(1- 6) 



n.2 



V22' 



^ "2 



^92n2(^0,0) 



-V, 



22 



— q'2n2(^0,0) 



+Op(l). 



(7) 



Substitute equations ([6]) and (171) into equation ([5]), we have 
2£(0,A)-2£(0o,Ao) 



n 



q2n2(^0,0) 



n 



q2n2(^0,0) 



T 



n 



\/n 



J 



J 



+ 



92n2(^0,0) 



T 



Hence the hmiting distribution of 2[l(6, \) — 1{6q,\q)] is the same as the 
distribution of ( [j^ ^ ^ f/^ ) under i/o- 

□ 
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Table 2: MSE Ratio of RSPELE to MLE with the smoothing parameter h = n2 to the 
power of -1, -3/4, -1/2 and -1/4. 



Distribution of sample 2 
n2iV(0, 1)^^(0,1-25) iV(0, 1-5) iV(0, 2) A^(0, 3) ^5 DE{0,0.5) DE{0,1) DE{0,1.5) 



h = 



10 


0, 


.739 


0, 


.841 


0, 


.922 


1, 


.019 


1 


.039 


0.747 


0. 


.719 


0, 


.364 


0, 


.664 


0. 


,847 


20 


0, 


.539 


0, 


.658 


0, 


.759 


0, 


.863 


0, 


.955 


0.609 


0. 


.560 


0, 


.187 


0, 


.487 


0. 


,697 


30 


0, 


.432 


0, 


.555 


0, 


.649 


0, 


.770 


0, 


.922 


0.442 


0. 


.437 


0, 


.120 


0, 


.362 


0. 


,572 



h — '^^^ 



10 


0, 


,700 


0, 


.810 


0, 


.878 


0, 


.987 


1 


.036 


0.730 


0, 


.692 


0, 


,339 


0, 


.645 


0, 


.821 


20 


0, 


,503 


0, 


.624 


0, 


.734 


0, 


.853 


0, 


.941 


0.588 


0, 


.512 


0, 


,172 


0, 


.466 


0, 


.662 


30 


0, 


,405 


0, 


.528 


0, 


.624 


0, 


.765 


0, 


.892 


0.418 


0, 


.413 


0, 


,107 


0, 


.333 


0, 


.552 



10 


0, 


.666 


0, 


.780 


0, 


.856 


0, 


,959 


1, 


.038 


0.702 


0, 


.665 


0, 


,326 


0, 


.625 


0. 


,795 


20 


0, 


.453 


0, 


.579 


0, 


.674 


0, 


,802 


0, 


.916 


0.543 


0, 


.475 


0, 


,162 


0, 


.433 


0. 


,616 


30 


0, 


.356 


0, 


.476 


0, 


.580 


0, 


,726 


0, 


.855 


0.385 


0, 


.369 


0, 


,101 


0, 


.319 


0, 


,524 



h = 712 



10 


0, 


.616 


0, 


.752 


0, 


.842 


0, 


.943 


1 


.036 


0.664 


0, 


.630 


0, 


,329 


0, 


.630 


0, 


,779 


20 


0, 


.403 


0, 


.532 


0, 


.635 


0, 


.773 


0, 


.899 


0.504 


0, 


.438 


0, 


,174 


0, 


.419 


0, 


,600 


30 


0, 


.320 


0, 


.434 


0, 


.532 


0, 


.682 


0, 


.838 


0.356 


0, 


.338 


0, 


,111 


0, 


.315 


0, 


,506 
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Table 3: Coverage Probability of MLE with the distribution of seeond population Normal, 
ts, ^5 and Double Exponential, numbers in the brackets are AL, a = nominal level. 

= U.8UU u.yuu u.ybu u.yyu u.8uu u.yuu u.ybu u.ygir 



n2 








Normal 








t3 


10 


0. 


,774 


0, 


.864 0.913 


0, 


.972 


0. 


,774 0.854 0.915 0.973 




(0. 


,774) 


(0, 


.994) (1.184) 


(1 


.557) 


(0. 


,777) (0.998) (1.189) (1.563) 


20 


0. 


,773 


0, 


.871 0.913 


0, 


.976 


0. 


,769 0.849 0.901 0.964 




(0. 


,784) 


(1. 


.007) (1.200) 


(1 


.577) 


(0. 


,777) (0.997) (1.188) (1.562) 


30 


0. 


,810 


0, 


.888 0.y32 


0, 


.983 


0. 


,760 0.858 0.912 0.972 




(0. 


,788) 




.011) (1.204) 


(1 


.583) 


(0. 


,786) (1.009) (1.203) (1.581) 










^5 








Double Exponential 


10 


0. 


,737 


0, 


.858 0.913 


0, 


.966 


0. 


,771 0.863 0.913 0.962 




(0. 


,789) 


(1. 


.013) (1.207) 


(1 


.586) 


(0. 


,781) (1.002) (1.194) (1.569) 


20 


0. 


,759 


0, 


.862 0.915 


0, 


.975 


0. 


,771 0.863 0.913 0.962 




(0. 


,795) 


(1. 


.020) (1.216) 


(1 


.598) 


(0. 


,781) (1.002) (1.194) (1.569) 


30 


0. 


,769 


0, 


.880 0.928 


0, 


.977 


0. 


,771 0.863 0.913 0.962 




(0. 


,790) 


(1. 


.015) (1.209) 


(1 


.589) 


(0. 


,781) (1.002) (1.194) (1.569) 
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Table 4: Coverage Probability of RSPELE with different distributions of second population 
by using median estimating equation, numbers in the brackets are AL, a ~ nominal level. 



a U.8UU u.yuu u.ybu u.yyu 



U.8UU u.yuu u.ybu u.yyu 



U.8UU u.yuu u.ybu u.yyu 



n2 






7V(0,1) 






10 





.754 0, 


846 902 





.963 




(o: 


.725) (0, 


.930) (1.102) 


(1 


.435) 


20 


0, 


.771 0, 


.868 0.913 





.964 




(0, 


.708) (0, 


.910) (1.085) 


(1 


.415) 


30 


0, 


.805 0, 


.888 0.935 





.976 




(0, 


.679) (0, 


.874) (1.047) 


(1 


.389) 








iV(0, 1.5) 






10 


0, 


.748 0, 


.845 0.899 





.962 




(0. 


.735) (0, 


.945)(1.116) 


(1 


.457) 


20 


0, 


.761 0, 


.864 0.911 





.961 




(0. 


,731) (0, 


.938)(1.111) 


(1 


.447) 


30 


0, 


.799 0, 


.882 0.936 





.973 




(0. 


.715)(0, 


.912) (1.084) 


(1 


.422) 








7V(0,2) 






10 


0, 


.748 0, 


.842 0.899 





.960 




(0. 


,741) (0, 


.948) (1.120) 


(1 


.464) 


20 


0, 


.753 0, 


.858 0.909 





.958 




(0. 


,743) (0, 


.951)(1.127) 


(1 


.462) 


30 


0, 


.798 0, 


.879 0.931 





.972 




(0. 


,736) (0, 


.936)(1.111) 


(1 


.456) 







iV(0,2) 




0. 


,746 0, 


.839 0.898 


.960 


(0. 


,744) (0, 


.949) (1.126) (1 


.471) 


0. 


,756 0, 


.855 0.908 


.958 


(0, 


,754) (0, 


.961)(1.138)(1 


.482) 


0. 


,791 0, 


.873 0.925 


.970 


(0. 


,754) (0, 


.957)(1.135)(1 


.483) 



0.758 0.848 0.908 0.964 
(0.733)(0.938)(1.112)(1.460) 

0.766 0.851 0.892 0.951 
(0.711)(0.911)(1.080)(1.411) 

0.771 0.869 0.916 0.966 
(0.700) (0.897) (1.068) (1.402) 

h 

0.734 0.840 0.897 0.963 
(0.744) (0.949) (1.124) (1.467) 

0.761 0.858 0.906 0.974 
(0.716) (0.916) (1.096) (1.441) 

0.766 0.877 0.929 0.971 
(0.693) (0.892) (1.065) (1.403) 





DE(0, 0.5) 




n 


759 854 909 


951 


(0, 


.716)(0.919)(1.094)(1 


.441) 


0, 


.780 0.872 0.918 


.966 


(0, 


.678)(0.878)(1.055)(1 


.401) 


0, 


.780 0.883 0.920 


.971 


(0, 


.643) (0.843) (1.018) (1 


.368) 








0, 


.758 0.848 0.907 


.951 


(0, 


.736)(0.941)(1.119)(1 


.460) 


0, 


.770 0.859 0.912 


.962 


(0, 


.711)(0.911)(1.086)(1 


.425) 


0, 


.772 0.866 0.915 


.965 


(0, 


.685) (0.882) (1.053) (1 


.392) 




DE{0,2) 




0, 


.755 0.847 0.904 


.951 


(0, 


.743)(0.952)(1.132)(1 


.470) 


0, 


.762 0.859 0.908 


.963 


(0, 


.731)(0.935)(1.109)(1 


.445) 


0, 


.770 0.866 0.906 


.961 


(0, 


.712)(0.913)(1.086)(1 


.416) 
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Table 5: Coverage Probability of RSPELE with different distributions of second popula- 
tion, and the smoothing parameter h = n.^^, numbers in the brackets are AL, a = nominal 

level. 

a 0.800 0.900 0.950 0.990 0.800 0.900 0.950 0.990 0.800 0.900 0.950 0.990 



n2 






7V(0,1) 










N(0,2) 






i:>-B(0, 0.5) 




1 n 


n 


7fi3 n 


861 920 


n 


973 


n 


752 


.840 0, 


.902 


966 


n 


787 n 889 934 


979 




(0, 


.687) (0, 


.886) (1.058) 


(1 


.394) 


(0, 


.754) (0, 


.963) (1, 


.139)(1 


.483) 


(0, 


.626) (0.829) (1.009) (1 


.368) 


20 


0, 


.797 0, 


.875 0.915 





.968 


0. 


.765 0, 


.863 0, 


.907 


.956 


0. 


.813 0.903 0.941 


.979 




(0, 


.649) (0, 


.840) (1.006) 


(1 


.344) 


(0, 


.756) (0, 


.963) (1, 


.138)(1 


.477) 


(0, 


.568) (0.766) (0.945) (1 


.311) 


30 


0, 


.816 0, 


.912 0.941 





.980 


0. 


.802 0, 


.885 0, 


.926 


.968 


0. 


.821 0.901 0.946 


.979 




(0, 


.610) (0, 


.796) (0.965) 


(1 


.305) 


(0, 


.751)(0, 


.956)(1, 


.130)(1 


.478) 


(0. 


.523) (0.718) (0.892) (1 


.263) 








iV(0, 1-5) 


















DE{0,1) 




10 


0, 


.756 0, 


.855 0.914 





.972 


0, 


.782 0, 


.873 0, 


.923 


.968 


0, 


.777 0.869 0.920 


.965 




(0. 


.726) (0, 


.930) (1.104) 


(1 


.448) 


(0, 


.709) (0, 


.910)(1, 


.082)(1 


.441) 


(0, 


.703) (0.909) (1.084) (1 


.438) 


20 


0, 


.786 0, 


.873 0.913 





.966 


0. 


.774 0, 


.869 0, 


.911 


.958 


0. 


.794 0.875 0.927 


.970 




(0, 


.707) (0, 


.908) (1.077) 


(1 


.415) 


(0, 


.668) (0, 


.861)(1, 


.028) (1 


.361) 


(0, 


.660) (0.856) (1.029) (1 


.370) 


30 


0, 


.803 0, 


.899 0.935 





.974 


0, 


.793 0, 


.880 0, 


.924 


.971 


0, 


.785 0.878 0.926 


.971 




(0. 


.684) (0, 


.874) (1.041) 


(1 


.377) 


(0, 


.648) (0, 


.836) (1, 


.001)(1 


.341) 


(0. 


.622)(0.811)(0.981)(1 


.320) 








7V(0,2) 










^5 








DE{0,2) 




10 


0, 


.752 0, 


.846 0.910 





.970 


0. 


.765 0, 


.865 0, 


.911 


.964 


0. 


.774 0.862 0.915 


.958 




(0. 


.743) (0, 


.948) (1.123) 


(1 


.473) 


(0. 


.716) (0 


.914) (1, 


.091)(1 


.436) 


(0. 


.732)(0.943)(1.121)(1 


.467) 


20 


0, 


.771 0, 


.871 0.909 





.961 


0. 


.785 0, 


.881 0, 


.922 


.976 


0. 


.780 0.868 0.919 


.966 




(0. 


.732) (0, 


.936)(1.112) 


(1 


.453) 


(0. 


.666) (0, 


.855)(1, 


.028)(1 


.376) 


(0. 


.707) (0.907) (1.078) (1 


.417) 


30 


0, 


.798 0, 


.893 0.934 





.973 


0, 


.799 0, 


.893 0, 


.939 


.974 


0, 


.780 0.874 0.917 


.965 




(0. 


.720) (0, 


.917) (1.087) 


(1 


.433) 


(0, 


.637) (0, 


.822)(0, 


.990) (1 


.334) 


(0, 


.677) (0.873) (1.043) (1 


.381) 
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Table 6: Coverage Probability of RSPELE with different distributions of second popu- 

— 1/2 

lation, and the smoothing parameter /i = rig , numbers in the brackets are AL. a = 

nominal level. 

a 0.800 O.yOU 0.950 0.990 0.800 0.900 0.950 0.990 0.800 0.900 0.950 0.990 



n2 




N{0,1) 








iV(0,2) 






DE(0,0.5) 




inn 


770 n 


.ooo u.yio u. 


.y / '4: 


u. 




.846 0, 


.908 0, 


070 


n 

u. 


iTX n 877 n Ok'vx n 


ORO 

.you 


(0, 


.630) (0 


.814) (0.977) (1, 


.310) 


(0, 


.756) (0, 


.965)(1, 


.143)(1, 


.491) 


(0, 


.508) (0.685) (0.861) (1, 


.240) 


20 0, 


.777 0, 


.870 0.921 0, 


.973 


0, 


.773 0, 


.865 0, 


.902 0, 


.954 


0, 


.797 0.899 0.945 0, 


.982 


(0, 


.551) (0, 


.713)(0.862)(1, 


.176) 


(0, 


.745) (0, 


.953)(1, 


.127)(1, 


.459) 


(0, 


.391)(0.531)(0.676)(1, 


.051) 


30 


.805 


.904 0.936 0, 


.980 


0, 


.800 0, 


.887 0, 


.933 


.972 


0, 


.795 0.893 0.943 0, 


.987 


(0, 


.491) (0, 


.636) (0.770) (1, 


.079) 


(0, 


.732) (0, 


,933) (1, 


.106)(1, 


.442) 


(0, 


.313) (0.431) (0.568) (0, 


.910) 






7V(0, 1.5) 








ti 












10 0, 


.763 


.849 0.911 0, 


.970 


0, 


.782 0, 


.874 0, 


.923 


.966 


0, 


.758 0.876 0.925 0, 


.969 


(0, 


.702) (0, 


.903) (1.074) (1, 


.416) 


(0, 


.662) (0, 


.853)(1, 


.021)(1, 


.375) 


(0, 


.649) (0.845) (1.019) (1, 


.377) 


20 0, 


.780 0, 


.864 0.911 0, 


.964 


0, 


.782 0, 


.862 0, 


.916 0, 


.967 


0, 


.791 0.889 0.923 0, 


.975 


(0, 


.657) (0, 


.844) (1.003) (1, 


.327) 


(0, 


.589) (0, 


.763) (0, 


.918)(1, 


.241) 


(0, 


.566) (0.740) (0.897) (1, 


.235) 


30 0, 


.816 0, 


.900 0.934 0, 


.976 


0, 


.797 0, 


.881 0, 


.928 0, 


.978 


0, 


.777 0.885 0.928 0, 


.978 


(0, 


.611)(0, 


.787) (0.939) (1, 


.258) 


(0, 


.538) (0, 


.698) (0, 


.845)(1, 


.152) 


(0, 


.493) (0.649) (0.795) (1, 


.123) 






iV(0,2) 








^5 












10 0, 


.753 0, 


.852 0.905 0, 


.970 


0, 


.766 0, 


.862 0, 


.920 0, 


.967 


0, 


.755 0.864 0.908 0, 


.967 


(0, 


.733) (0, 


.940)(1.113)(1, 


.461) 


(0, 


.662) (0, 


.854) (1, 


.021)(1, 


.365) 


(0, 


.708) (0.910) (1.087) (1, 


.446) 


20 0, 


.774 0, 


.865 0.906 0, 


.959 


0, 


.786 0, 


.880 0, 


.936 


.972 


0, 


.786 0.866 0.922 0, 


.963 


(0, 


.706) (0, 


.906) (1.073) (1, 


.402) 


(0, 


.578) (0, 


.747) (0, 


.900) (1, 


.226) 


(0, 


.649) (0.836) (1.002) (1, 


.342) 


30 0, 


.809 0, 


.899 0.935 0, 


.973 


0, 


.803 0, 


.901 0, 


.940 0, 


.975 


0, 


.785 0.875 0.919 0, 


.968 


(0, 


.677) (0, 


.862) (1.026) (1, 


.359) 


(0, 


.520) (0, 


,679) (0, 


.823)(1, 


.137) 


(0, 


.596) (0.774) (0.934) (1, 


.264) 
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